Quantum statistical effects in nano-oscillator arrays 
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We have theoretically predicted the density of states(DOS), the low temperature specific heat, 
and Brillouin scattering spectra of a large, free standing array of coupled nano-oscillators. We have 
found significant gaps in the DOS of 2D elastic systems, and predict the average DOS to be nearly 
independent of frequency over a broad band / < 50GHz. At low temperatures, the measurements 
probe the quantum statistics obeyed by rigid body modes of the array and, thus, could be used to 
verify the quantization of the associated energy levels. These states, in turn, involve center-of mass 
motion of large numbers of atoms, N > 10 14 , and therefore such observations would extend the 
domain in which quantum mechanics has been experimentally tested. We have found the required 
measurement capability to carry out this investigation to be within reach of current technology. 



Over the past decade, significant effort has been fo- 
cused on the observation of quantum effects in large sys- 
tems, both to test quantum mechanics in this regime 
and to directly observe quantum decoherence phenom- 
ena arising from coupling to the environment. Mechan- 
ical oscillators have often been considered in this con- 
text because displacement states of the oscillator can 
serve as "Schrodinger cat" states in the investigation of 
these phenomena. The measurement concepts reported 
to date in this areapj S LH 0] have primarily involved 
investigating the quantum properties of a single oscil- 
lator, though short linear chains of coupled oscillators 
have been considered recently in the context of minimum 
thermal conductance and entanglement Q. Concepts 
involving a mechanical oscillator coupled to the optical 
modes of a cavity have been considered by Bose et. al.Q, 
while the investigation of a mechanical oscillator coupled 
to the electrons in a Cooper pair box 9] , the electronic 
analogy, is underway. Recent experimental results have 
pushed close to the limits of the uncertainty principle (lcj. 

A different approach, and the central focus of this pa- 
per, is to explore the extent to which one can probe the 
quantum behavior of mechanical oscillators by observ- 
ing the effectively bulk properties of a large, free stand- 
ing array of coupled nano-oscillators. Measurements of 
this type avoid the significant experimental difficulties 
that direct observations of the motion of individual nano- 
mechanical oscillators involve 0, 0, ^| and offer a new 
avenue to explore their quantum behavior. The simplest 
bulk measurements do not explore the coherent quantum 
behavior of the system, but they do enable us to ver- 
ify that macroscopic whole body modes obey quantum 
statistics and therefore have quantized energy levels. 

The quantum behavior of the system is expected to be 
governed by the field operators 
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are the classical eigenfunctions. We have verified this in- 
tuitive result via the canonical quantization prescription, 
but the details are too lengthy to be given here; for our 
purposes in this letter, we assume Eq. 0} to be valid. 

The classical solution for the eigenfunctions and the 
spectrum is evidently central to the quantum behavior 
of the system (e.g. the presence and size of frequency 
gaps and the nature of the eigenfunctions). Typically, 
the dynamic behavior of even a simple elastic structure 
contains significant complexity beyond a model of cou- 
pled simple oscillators [llj. Therefore, below, we give nu- 
merical results for the classical modes of a generic, peri- 
odic, nano-structured system. We have found large gaps 
in the spectrum over a broad frequency range, .1GHz 
< f < 50GHz. We then obtain theoretical predictions 
for two of the most basic bulk measurements of such a 
system, the specific heat and Brillouin scattering. At 
fairly low temperatures, T ~ 30mK, the quantum statis- 
tical behavior of modes involving whole body motion of 
the nano-oscillator array is observable. 

The geometry we consider is that of a periodic, planar 
structure as shown in Fig. 1. The particular dimensions 
chosen can, of course, be rescaled by a factor A accord- 



ing to L — > XL, 



We assume, for simplicity, 




where the 0^(00) are the usual creation (annihilation) op- 
erators for mode a, p is the mass density, and the ip a (x) 



FIG. 1: Geometry of a periodic planar nano-structured 
phononic crystal. The values used for calculations are: / = 
200nm, a = 360nm, h — 20nm, thickness = 20nm, and mate- 
rial parameters of Si. 
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that the material is isotropic. The normal modes are cal- 
culated by seeking solutions of the Bloch form ip(x) = 
exp(ikx)U (x) , where k = (k x ,k y ) = (2mr/L x , 2mir / L y ) 1 
and U (x) is assumed to possess the periodicity a of the 
lattice. Here, n and m are integers, L x and L y are the 
dimensions of the crystal, and the boundary conditions 
are taken to be periodic, an approximation which will 
not be significant for a large system. This results in a 
boundary value problem in the unit cell which yields the 
resonance frequencies for a particular wavenumber, with 
one frequency for each sub-band of the spectrum. 

The numerically computed density of states(DOS) is 
shown in Fig. 2. The vibrational modes that ex- 
ecute whole body motion locally are the modes of 
greatest interest with regard to testing the limits of 
quantum mechanics^^ an d environmentally induced 
decoherence0,0|, because the associated energy eigen- 
states exhibit coherent center of mass motion of large 
numbers of atoms. The states in the lowest frequency 
band extending up to about 2GHz are in large part whole 
body modes involving significant domains of the entire 
structure. The first quasi-gap in the spectrum, associ- 
ated with weakly coupled local 1-1 modes (the peak at « 
2.25GHz), occurs in the range from l-2GHz, while a true 
gap appears in the spectrum at approximately 6GHz. 

The presence of gaps is significant because modes near 
the band edges are Anderson localized by irregularity 
This phenomena may be advantageous in the experimen- 
tal study of the quantum behavior of the system be- 
cause of the significantly larger physical displacements 
for the same energy levels. Further, localized wave func- 
tions vanish exponentially on the boundaries and there- 
fore promise to suffer much less decoherence resulting 
from coupling to the environment. 

A notable feature of the results is that the average DOS 
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FIG. 2: Density of phonon states of a periodic planar nano- 
structured array. Solid: Total DOS; dashed: average DOS. 
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FIG. 3: Specific heat of a periodic planar nano-structured 
array. 

is more or less constant in this frequency range. This be- 
havior results from the dispersive properties of the flexu- 
ral modes which dominate the spectrum because of their 
slow phase velocity. One may see this qualitatively by 
estimating the DOS of a uniform, flexural wave system. 

The dispersion relation for flexural waves propagating 
in a thin plate of thickness t is ui = (D/m) 1 ^ 2 k'j, with 
D = Et 3 /[12(1 - v 2 )\ the flexural rigidity, and m = pt 
the areal density. Here E is the Young's modulus, p the 
density, c/ = w/kf the flexural wave speed, and v the 
Poisson ratio. The DOS of a uniform, d-dimensional, 
flexural wave system of size L d is therefore given by 

and is independent of frequency in two dimensions. 

Theoretical predictions for the specific heat can be ob- 
tained directly from the DOS/vol, n(uS), via 

du k b f e*M 

Cv = — — I dui n(uo)- — t- x — [A) 

dT V J v '(e*M - l) 2 w 

where x(u>) — hwjk^T. The result is shown in the tem- 
perature range below 250mK in Fig. 3a. In this range, 
the specific heat depends on the temperature in an ap- 
proximately linear fashion. This behavior differs from the 
typically quoted T 2 dependence 0] of the specific heat 
in 2-D phonon systems because of the roughly constant 
DOS at low frequencies, an aspect discussed above. 

Employing Eq. (J2J , one may show that the low temper- 
ature heat capacity of a uniform d-dimensional flexural 
system of size L d is 

where £ is the Riemann zeta function. In 1-D, the heat 
capacity varies as T 1 / 2 , while in 2-D the predicted tem- 
perature dependence is T, as we have found above. The 
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heat capacity of low dimensional phonon systems is there- 
fore larger at low temperatures than predicted using the 
customary T 2 temperature dependence of C v . 

The nano-structure of the system, as opposed to the 
overall dimensionality, has a relatively weak impact on 
the temperature dependence of the specific heat, giving 
rise to the bump in the vicinity of T ~ 35mK . Neverthe- 
less, the observation of a linear temperature dependence 
at low temperatures would provide direct evidence that 
the relevant whole body modes are behaving as simple 
quantum oscillators obeying Bose statistics. 

Experimentally, a 1000 x 1000 array of oscillators with 
dimensions as given in Fig.l would result in a heat ca- 
pacity of 1.6 x 10~ 16 J/K at 100 mK, an estimate which 
does not include the addenda (heater and thermometer) 
necessary for the measurement. This is a challenge to 
the current limit of 2 x 10- 9 J/K at 1.5K [H0- But 
with the rapid advance of nano-technology and the lower 
temperature range providing more sensitivity, it is not 
unachievable. 

The quantum statistics obeyed by the whole body 
modes of a nano-oscillator array can also be probed in 
an optical scattering experiment, and the results so ob- 
tained would be preferable to those from a specific heat 
measurement, enabling one to observe the behavior of 
particular modes. We approximate the scattered field 
due to thermal vibration of a nano-structured array by 
estimating the classical electromagnetic scattering from a 
vibrating structure, and then replacing the required aver- 
ages over displacement fields in the classical expressions 
by quantum expectation values. The scattered compo- 
nent of the vector potential is given by 



A{x,u) = ^ I J eq {x") 
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where the equivalent source current is J eq — — iw(e — 
eo)E= —iwtoxE with e the dielectric constant, x the sus- 
ceptability of the system, and E the total electric field 
inside the the dielectric. To leading order, E can be re- 
placed by the incident field in the integration over the 
source, a Born approximation which is reasonable be- 
cause of the small thickness of the structure. 

The evaluation of the scattering cross section can be 
further simplified by the following approximations. We 
assume: the motion is slow in comparison with optical 
frequencies; the displacements and sample thickness are 
small relative to an optical wavelength; and, finally, the 
displacements are dominated by the normal components. 
One may then show that the cross section is given by 
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Here q = k — ki is the change of the optical wavenumbcr 
and Q is the change of the incident frequency. The quan- 
tity ao(<?) = J d 2 xexp(—iqx) is a static aperture integral 



over the surface and gives the Rayleigh scattering, while 
ip(q, fl) is the normal displacement field of the system, 
and gives the dynamic response. A reasonable approxi- 
mation to the cross section is then obtained by replacing 
\ip(q, ri)| 2 by the thermal average of the quantum expec- 
tation value of Eq. 1, 



<\^(q,n)\ 2 > T = [da] 
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where n(uj ai T) is the Bose distribution, and the two 
terms correspond to Stokes and anti-Stokes peaks respec- 
tively in the usual way for Brillouin scattering. For ex- 
tended modes, the wave function tp a takes on the sim- 
ple form, ip n k{q) = (2n/a) 2 5 2 N (q - k)l 2 j nk (q - k), where 
5'jq{k) = exp(— zfc-ma) is a lattice delta function, and 
j n k{q) = l~ 2 J d 2 x exp(— iqx)u n k(x) is the modal optical 
source strength, the normalized Fourier transform of the 
mode shape of a single element of the array. 

Consider the spectra associated with a phonon of 
wavenumber q which scatters light at angle sydT 1 {q / 2k) . 
Using Eq. JHJ and Eq. Q, we estimate the ratio R 
of the Brillouin peak to Rayleigh scattered levels for a 
100/imx 100/im sample. For normal incidence, scattering 
angle 9, and 1GHz phonons, we find R ~ 4 x 10~ 6 # 4 
away from the main Rayleigh lobe in directions along a 
diagonal (9 = <fi = 7r/4). Using a multi-pass Fabry-Perot 
interferometer with contrast ratio 10 9 , the Brillouin peak 
level is found to equal the background (and thus be read- 
ily observable) for 9 > 7 deg. 

Having established that the Brillouin spectra is ob- 
servable above the Rayleigh background, we consider the 
scattered power level. The total power scattered into 
solid angle AQ near an anti-Stokes peak may be obtained 
from Eqs. © and (JjJ as 
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n(u) q ,T)An . (8) 



Here Iq is the total power incident on the array; k is the 
optical wavenumber; r — (kdx/2) is the plane wave re- 
flection coefficient; j„ 9 (0) is the optical source strength; 
Ax — (h/2pl 2 uj q ) 1 / 2 is the size of the ground state 
wave packet of a single element of the array; and Af2 
is the solid angle in which light is gathered. The de- 
pendence of the scattered power on the extent of the 
ground state wave function for a single clement of the ar- 
ray is a result of coherent scattering associated with the 
N oscillators, each moving with the average displacement 
Axm = (ti/2Npl 2 uj q ) 1 / 2 . The peak scattered Brillouin 
power at low temperatures therefore scales as N. 

Note also that the scattered power is proportional 
to the mean number n(uj qi T) (or (n(u> q ,T) + 1) for 
the Stokes peak) of phonons present at temperature T. 
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Therefore, a measurement of the temperature depen- 
dence of the Stokes (anti-Stokes) peak would provide di- 
rect evidence that the mode q is obeying quantum statis- 
tics. This requires measurements in the temperature 
range kT ~ Tioj q . 

We now give an order of magnitude estimate of the 
power scattered into a Stokes peak for the geometry pre- 
viously considered at low temperatures. Using the same 
parameters as employed above and, further, assuming 
ink = AS1 = 1, we find / ~ 4 x 10~ 14 /o at low temper- 
atures such that n{bJ q ,T) = 1. Taking I — 100/zW, a 
heat load well within the cooling capacity of many cur- 
rent systems, we obtain a Brillouin scattered power of 
approximately 10 photons/sec. This is indeed a small 
signal, but note that the power level I in Eq. Q scales 
as A 2 if we rescale the nano-structure, L — » XL, provided 
we do not violate the small dimension approximations. In 
obtaining this scaling law, we have assumed the incident 
power to be fixed; if, instead, the power flux is regarded 
as fixed, an additional factor of A 2 results. Thus, higher 
signal levels will be observed for larger structures, at the 
cost of lower temperatures. 

Topological disorder in the structure reduces the coher- 
ence of the Bloch modes and thus the levels of the Stoke's 
and anti-Stoke's peaks. Near band edges, however, the 
Bloch modes of the array are localized by irregularity in 
the system^5|, and the higher resulting vibration lev- 
els partially compensate for the loss of coherence. The 
mean square displacement level associated with quantum 
level motion increases as £~ 2 as the localization length 
£ decreases. Assuming the q = 2\k\ sin 0/2 selection rule 
breaks down as a result of the localization, as in the case 
of amorphous materials^^ , all q's contribute to the scat- 
tering at each angle, and the strength of the scattering 
is proportional to the area of mode coherence I 2 (I is the 
coherence length). The cross section for Brillouin scatter- 
ing will, therefore, be proportional to {l/^) 2 n{uo), where 
n(u) is the density of states. For strong localization of a 
resonance near a band edge, l/£ is of order 1, and n(u>) 
is of order N, leading to a scattered power of the same 
order of magnitude as given in Eq. (JHJ . This phenomena 
should be observable near the band edges associated with 
the gaps at 1.5GHz and 6GHz for the geometry consid- 
ered here. 

We have explored the extent to which the quantum 
behavior of mechanical oscillators can be probed by ob- 
serving two of the simplest physical properties of a large, 
free standing array of coupled nano-oscillators: specific 
heat and Brillouin scattering. In both cases, the mea- 
surements are challenging but within reach of current 
capability. Carrying out such experiments at low temper- 
atures would enable us to observe the quantum statistics 



obeyed by rigid body modes involving large numbers of 
atoms and, thus, verify that the energy levels of these 
states are quantized. Such measurements would signifi- 
cantly enlarge the domain in which quantum mechanical 
predictions have been tested experimentally. We have 
not examined the possibility of employing such a system 
to explore coherent quantum phenomena, but because 
the observation of small occupation number states is not 
out of reach, such an attempt seems promising. 

The authors gratefully acknowledge the Office of Naval 
Research for support of this work. 
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